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We analyze critical behavior of superfluid density Ps in strongly disordered superconductors near 
superconductor-insulator transition and compare it with the behavior of the spectral gap A for col¬ 
lective excitations. We show that in contrast to conventional superconductors, the superconductors 
with preformed pairs display unusual scaling relation pa oc close to superconductor-insulator 
transition. This relation have been reported in very recent experiments. 


In a disordered conductor the superconductivity is first 
suppressed and then completely destroyed at high dis¬ 
order. There are two mechanisms for the suppression 
(and eventually the full destruction) of the superconduc¬ 
tivity by disorder (for recent reviews se^ ancPl, chapter 
1). The first mechanism attributes the suppression of 
the superconductivity to the increase of the Coulomb in¬ 
teraction that results in the decrease of the attraction 
between electrons and their eventual depairing.l^ In this 
mechanism the state formed upon the destruction of the 
superconductor is essentially a poor conductor. The al¬ 
ternative mechanism attributes superconductivity sup¬ 
pression to the localization of Cooper pairs that remain 
intact even when superconductivity is completely sup¬ 
pressed. The latter mechanism is called bosonic and the 
former fermionic. The theory of the bosonic mechanism 
has a long history: this scenario of the superconductor- 
insulator transition was suggested long agcPUl but was 
not developed further until recentljP^ when experimen¬ 
tal datcPSS indicated the existence of a few materials 
that show such behavior. The bosonic mechanism is also 
supported by numerical computation^^. 

The interest to the superconductor-insulator transition 
without Cooper pair destruction is both fundamental and 
practical. First, it provides a perfect example of the dis¬ 
order driven quantum transition in the closed system. 
Second, the bosonic superconductor in the vicinity of the 
transition might be ideal element for the isolation of the 
coherent quantum system from the environmenl!i^ and 
for sensitive detectors of microwave radiatiorP^. One of 
the most important properties for these applications is 
the value of the superfluid density of the superconduc¬ 
tor. The goal of this paper is to compute this quantity 
in the bosonic mechanism; we also show that its behav¬ 
ior as the superconductivity is suppressed distinguishes 
fermionic and bosonic mechanisms of the superconduc¬ 
tivity suppression. 

Characteristic feature of the fermionic mechanism is 
that both the transition temperature Tc and the spec¬ 
tral gap A at T <C Tc are suppressed simultaneously, so 
that their ratio is left nearly constant in the broad range 
of Tc variation. At all disorders the superconductor re¬ 
mains qualitatively similar to a conventional BCS super¬ 
conductor. When the superconductivity is completely 
suppressed by disorder the normal metal is formed, per¬ 


haps with a weak tendency towards localization. 

In contrast the state formed in the bosonic mecha¬ 
nism is qualitatively different from both the normal metal 
and conventional superconductor. In this mechanism 
the superconductivity disappears at the Superconductor- 
Insulator transition (SIT) whilst the gap in the single 
electron spectrum remains intact!^. Superconducting 
state formed in the vicinity of SIT is therefore a pseu¬ 
dogapped stat^. The best example of such behavior is 
amorphous InO films in which the critical temperature 
decreases by more than a factor of three while the single¬ 
particle gap stays practically independent on Tc but fluc¬ 
tuates spatially (Ai « 0.45 — 0.6 meV)!^. Another evi¬ 
dence of the bosonic mechanism is formation of a strong 
insulator characterized by a near ly-act ivated dependence 
of resistivity R{T) as observed irPSEH. 

The bosonic mechanism attributes the suppression of 
the superconductivity to the competition between Ander¬ 
son localization and Cooper attraction between electrons 
assuming that the Coulomb interaction plays a minor 
role. The key feature of this mechanism is that the su¬ 
perconducting state is formed by the electrons in the lo¬ 
calized single-electron eigenstates ifiir), with a relatively 
large localization length Lioc which depends on the prox¬ 
imity of the Fermi-energy to the Anderson mobility edge 
Ec- The presence of a length scale L/oc that is not related 
to superconducting properties of the conductor leads to 
the appearance of new energy scale in the problem that 
differs from the single electron gap in a usual BCS su¬ 
perconductor. Namely, attraction between two electrons 
that populate the same orbital state leads to the 

formation of a bound pair with a binding energy 

2Ai = dVV-f(r) ~ gL"? 

where g is the Cooper attraction constant (with dimen¬ 
sionality [Energy] X [Volume]), and ^2 « 1.3 is the fractal 
dimension of 3D Anderson transition, see 

The energy scale Ai discriminates against the odd (sin¬ 
gle electron population) while even (zero or two electrons) 
parts of the Hilbert space remain at low energies. As a 
result, single-particle density of states (DoS) acquires a 
pseudogap irrespectively of the presence of superconduct¬ 
ing correlations. Experimentally, the appearance of the 
pseudogap in the absence of long range superconducting 
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Figure 1. Sketch of the physical properties of the supercon¬ 
ductors close to the disorder driven superconductor-insulator 
transition. As the disorder is increased the transition tem¬ 
perature drops while single particle gap remains constant re¬ 
sulting in the deviations from BCS relation 2/S.iIT mf ~ 3.5 
(indicated by the dashed line). At the same time the super- 
fluid density p(0) at zero temperature also drops with the 
respect to the predictions of the BCS theory as shown by thin 
(red) line. The convenient measure of the disorder is pro¬ 
vided by the minimal conductance, Omin above the transition 
temperature. To make the connection with the experimental 
situation we have included in this plot the actual experimen¬ 
tal data froirP^ for InO. For these points a min is measured 
in mn~^cm~^. Very close to superconductor-insulator tran¬ 
sition superfluid density might be controlled by non-trivial 
critical exponents. 


order was demonstrated in InCP^ where strong suppres¬ 
sion of low-voltage tunneling conductance was found up 
to T > 2Tc. Similar data were reported for sufficiently 
disordered TiN filmJi^. In a pseudogapped supercon¬ 
ductor long-range correlations develop due to coupling 
between localized bound pairs, and lead to formation of 
another energy scale, A, that is related to superconduct¬ 
ing long range order. The appearance of this energy scale 
was observed re cently by Andreev tunneling spectroscopy 
in InOa; sampleJi^^^close to SIT with Tc l.hK . These 
data indicate the appearance of the the collective energy 
scale A < Ai that obeys BCS-like temperature depen¬ 
dence while Ai is almost temperature independent at 
T < Tc- Another demonstration of the presence of some 
smaller energy scale that lies within the single-particle 
gap is provided by optical spectroscopy of InOa, and NbN 
samples of various disorder.^ 

Thin superconducting films are known to possess one 
more important energy scale, 0, that characterizes the 
dependence of the free energy on the phase gradient: 

The helicity modulus 0 determines the location of the 
Berezinsky-Kosterlitz-Thouless (BKT) transition at 
the transition Q{Tbkt) = ^Tbkt- Energy 0 is related 
to the superfluid density Ps which determines supercur¬ 
rent in the film of thickness d in presence of a vector 


potential: j = —psAjc, namely 0 = (?i/2e)^ps • d. 

For conventional disordered superconductors there is a 
linear relation (at T T^) between superfluid density p^ 
and the energy gap A: 


Ps 


ttctA 

h 


( 1 ) 


where a is the normal-state resistivitji^. Although usu¬ 
ally derived in the framework of the BCS theory this 
relation is more robust because it can be traced to the 
optical weight conservation. So, not surprisingly, it 
is valid also for the moderately disordered InOx films 
with Tc ~ 2.5 — S.5K, as was found via the kinetic in¬ 
ductance measurement^. Magnetoresistance studies of 
these films, see 1^, as well as Andreev tunneling spec- 
troscopjl^, show expected properties of disordered BCS 
superconductors, namely: no pseudogap, Andreev gap 
that has the same value as the single-particle gap, fi¬ 
nally, the destruction of superconductivity by magnetic 
field leads to a normal metal state without noticeable 
R{B) peak. 

We now demonstrate that a pseudogapped supercon¬ 
ducting state is characterized by a different relation be¬ 
tween 0 and A, leading to a number of unique features. 
We start with the expression for the full spectral weight 
Rtot frequency-dependent conductivity as derived in 
Sec. 6.4 of Ref.® 





‘i’j 


AiAj tanh pSi tanh Psj 

Si€j 


( 2 ) 

Here P = I/T, V is the system’s volume, i,j enu¬ 
merate single-electron eigenfunctions those eigenvalues 
are matrix elements = J d^r'(/)|(r)'0|(r), and 

£i = \JQ + Af. The quantities Ai are the order pa¬ 
rameter amplitudes related to the superconducting order 
parameter in a real space A(r) by 


A(r) 


9 
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tanh(^^2(j.) 


( 3 ) 


In the mean-field approximation , the amplitudes A^ obey 
self-consistency equations 


A,= 


9 
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tanh(/3ej) 


( 4 ) 


Mean-field approximation of Ref.® assumes that ampli¬ 
tudes Ai are actually a slow functions of the single¬ 
particle energies ^i, i.e. one can replace Ai by a regular 
function A(^) evaluated at ^ this function has to 

be determined from the continuous version of Eq. 

A(e) = ^l dCM(C - (5) 


where A = gvo and vq is the DoS value in the normal 
state (per single spin projection), and M{uj) = VMij 
with - Cil = w- 
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Eq. (§ was derived under the assumption that eigen¬ 
functions ^ij (r) that contribute mostly to the sum over 
i,j are relatively well-localized, with the typical distances 
Xij ^ i?o between the maxima of their envelopes that are 
somewhat larger than localization length Lioc- We esti¬ 
mate effective interaction range as Rq ~ L/oc In ^, where 
5l = {’^oLfoc)~^ is the level spacing within localization 
volume. Pseudogaped superconducting state is realized 
when A < Ai ^ Sl, so that Rq ^ Ejoc- 

At low temperatures the optical sum rule arguments 
show that the major contribution to the total spectral 
weight comes from the superconducting density, i.e 
« /Os(0) = Ps-I^The Eq.([^ can be simplified by 
eliminating the sum over j with the help of Q and by 
replacing square of dipole matrix elements by its average, 
x^j Rq/2. At T = 0 we get 


Ps 


e^Rl V Af 

I 


2voe^Rl 2 


Jo x/F+W) 

( 7 ) 


We emphasize that A in this equation represents col¬ 
lective gap as measured by the Andreev spectroscopj^^ 
and in the THz optical measurement^^, which is smaller 
than single-particle gap Ai. 

The major difference between relations 0 and Q 
is that in the pseudogapped superconductor ~ 
whereas in the usual case ps ^ A. Note that general 
arguments related to optical weight conservation are not 
applicable for the pseudogapped superconductor due to 
the presence of the second energy scale, Ai. The unusual 
scaling of ps oc A^ in Eq.Q is due to the presence of an 
independent spatial scale Rq that determines the range 
of the tunneling matrix elements Mij between localized 
eigenstates K is crucial that Rq weakly depends 

on A. The counterpart of Rq in a usual dirty supercon¬ 
ductor is given by its low-temperature coherence length 
^0 ~ y^hD/A. In this case instead of R^ one should use 
(q oc D/a and linear relation pg oc A is restored. Note 
that quadratic scaling p^ oc A^ is known for superfluidity 
in Bose systems, see for exampl^^. 

The prediction of unusual scaling between A and ps is 
important close to SIT where both A and ps are expected 
to decrease strongly. Observation of the scaling ps oc 
A^ would serve as an independent demonstration of the 
bosonic nature of a superconductive state. 

The crucial step in the derivation of the result Q is 
the averaging over Ci and A^. Close to the transition to 
the insulating state the superconducting order parameter 
becomes very inhomogeneou^^, which requires reexam¬ 
ination of this procedure. In this regime distribution of 
the order parameter becomes very broad. In the Bethe 
lattice approximation it is given by 


R(A) 


Am 

^0 


^1-l-m 


Ao < A < A„ 


( 8 ) 


Here Aq denotes the typical (most probable) value of the 
local order parameter A, the exponent m is slightly less 


than unity: I —m = eA <C I where e = 2.718.... Distribu¬ 
tion (|^ is applicable up to the upper cut-off Amax Ai. 

Another kind of the order parameter distribution was 
obtained numerically by Lemarie et al for 2D attractive 
Hubbard model with strong random potentiaP^, it was 
claimed to be related to a distribution of the Tracy- 
Widom typ^^ that do not have power law tails of the 
form Q, but does have a large dispersion. 

In presence of strong statistical fluctuations of A^ the 
MEA equations (^1 are not valid, thus it is not possible 
to follow the route of calculations that lead to Eq.Q. 
Instead we note that original Eq. 0 contains double sum 
over sites i,j over large number ^ ^ 1 of statisti¬ 

cally independent terms for each i. Thus it is possible to 
estimate this sum by 


Ps 


voe^Rl A 2 


(9) 


where Aq is the typical (most probable) value of the local 
order parameter A^, irrespectively of its specific distribu¬ 
tion P{A). Qualitatively, one expects that rare pairs of 
states i,j with anomalously large A^ ^ Aq cannot con¬ 
tribute considerably to macroscopic superfluid density ps 
since they will be ’’screened” by weaker typical pairs due 
to redistribution of the supercurrent density. 

Notice that all quantities entering Eq. 0 except for 
the interaction range i?o are measurable. The latter can 
be determined if the scaling relation 0 is observed in a 
broad range of Aq and ps- 

Prediction 0 is in a rough qualitative agreement the 
data that became available very recentljED. Indeed, total 
evolution of the gap as measured in this experiments by 
optical spectroscopy, Eig.2c, spans about one order of 
magnitude between the ’’crossing point” corresponding 
to the reduced transition temperature Tj, = 0.5, and the 
most disordered samples with T). « 0.2. At the same 
time, superfluid density ps (extracted from the imaginary 
part of the optical impedance) changes by nearly 2 orders 
of magnitude in the same range of Tc, according to Fig.3b 
of the same paper. 

In the close vicinity of the quantum phase transition 
from superconducting to insulating state one expects crit¬ 
ical behavior characterized by the exponents that might 
differ from the mean field result (|^ , especially in low di¬ 
mension^ systems. For instance, the recent numerical 
workPflHU on the hard core boson model reported behav¬ 
ior Ps oc A“ with the exponent a « 2.5 in the critical 
regime. However, applicability of 2D scaling for strongly 
disordered superconducting films studied experimentally 
is not evident, because many of these films are not thin 
enough to be considered two-dimensional. 

Close to the critical point one expects large spatial fl- 
cutuations of the order parameter. The resulting fluctua¬ 
tions of the superfluid density were studied numerically in 
the recent papeJ^for the 2D hard-core boson model with 
random potential, using Quantum Monte Carlo method 
on systems with linear sizes in the interval L = 12 — 32. 
A broad probability distribution Vi/npa) was found both 
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in superfluid and insulating phases; however, on super¬ 
fluid side of the transition, the width of this distribu¬ 
tion diminishes with L, whereas an opposite tendency is 
found for the insulating state. These results indicate that 
in a macroscopic system dc superfluid density is a self¬ 
averaging quantity in the superconducting state. Analyt¬ 
ical calculation of the corresponding distribution 7^(In p^) 
in the framework of the theory developed here is an in¬ 
teresting problem which we leave for future studies. 

The spatial fluctuations of ps are probably observable 
in the low frequency measurements. Usually Ps is mea- 
su red v ia ac kinetic inductance C oc 1/ps, see for exam- 
pl (j32 | 33 | Nonzero measurement frequency w determines 
the length scale loj = \J~D]Z} where coherent transport 
takes place. With a typical diffusion contant for very 
strongy disordered superconductors D ~ 0.1 — lew? js, 
the length ranges from 0.02 — 0.5 micron for high- 
frequency measurment^^ to 10 — 100 micron for low- 


frequency one^^. In the latter case it might be possi¬ 
ble to detect spatial fluctuations of ps with a scanning 
technique. 

Sharp drop near SIT of the helicity modulus 0 = 
{T?/'le^)ps ■ d according to ^ leads to an enhancement 
of both thermal and quantum phase fluctuations, and to 
suppression of both and critical magnetic held with 
respect to the mean-field estimates, see also Ref.l^. 

In conclusion we have shown that the bosonic mecha¬ 
nism of the superconductor-insulator transition implies a 
nearly quadratic scaling of the superfluid density with the 
superconducting order parameter, the prediction that can 
be verified experimentally. This resulting very low values 
of the superfluid density would make materials close to 
the transition extremely useful for the applications. 
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